Photon production from an anisotropic quark-gluon plasma 
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We calculate photon production from a quark-gluon plasma which is anisotropic in momentum 
space including the Compton scattering, q{q) g — > g(g) 7, and annihilation, qq ^ g^, processes. We 
show that for a quark-gluon plasma which has an oblate momentum-space anisotropy, < (pr), 
the photon production rate has an angular dependence which is peaked transverse to the beam 
line. We propose to use the angular dependence of high-energy medium photon production to 
experimentally determine the degree of momentum-space isotropy of a quark-gluon plasma produced 
in relativistic heavy-ion collisions. 
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The ongoing heavy-ion experiments at the Relativistic 
Heavy Ion Collider (RHIC) and planned at the Large 
Hadron Collider (LHC) aim to determine information 
about QCD under extreme conditions. The goal of these 
experiments is to produce and study the properties of 
a quark-gluon plasma which is expected to be formed 
when the temperature of nuclear matter is raised above 
its critical value, Tc ~ 200 MeV. In the final state these 
experiments will generate on the order of 10"* hadronic 
particles. However, during the early stages of the sys- 
tem's evolution the produced "fireball" has an energy 
density on the order of 20 — 40 GeV/fm"^ and one should 
describe the system in terms of partonic (quasiparticle) 
degrees of freedom, quarks and gluons. The question of 
whether or not the partonic degrees of freedom interact 
frequently enough to be described as an isotropic ther- 
malized state before hadronizing and finally flying apart 
as non-interacting hadrons is an important and difficult 
question [1] with recent developments focusing on the 
impact of momentum-space anisotropics on plasma evo- 
lution [2-15]. 

Absent a theoretical proof that the system will be- 
come isotropic in momentum space (and possibly ther- 
malized) quickly enough, one is forced to look for observ- 
ables which would be sensitive to early time momentum- 
space anisotropies in the quark and gluon distribution 
functions. One such observable is high-energy medium 
photon production. Because photons interact weakly 
with the surrounding matter during the evolution of the 
plasma and subsequent hadronic phase they are ideally 
suited for determining information about the early stages 
of a heavy-ion collision. For this reason thermal pho- 
ton production from an isotropic quark-gluon plasma has 
been studied extensively [16-23]. Thermal photons make 
significant contributions to the total photon rate in the 
momentum range 2 GeV < Pt < 6 GeV. Here we are 



also interested in the medium produced photons in this 
momentum range; however, we relax the assumption of 
isotropy in order to see what effect this has on the angu- 
lar dependence of medium photon production. 

For our final results we present the energy and an- 
gular dependence of the photon production for two dif- 
ferent system anisotropies. The results show that for 
finite anisotropies photon production is peaked trans- 
verse to the anisotropy direction. For oblate geome- 
tries, (pj^) < (j)'ip), this corresponds to production being 
peaked transverse to the beamline. This result suggests 
that it may be possible to determine the time-dependent 
anisotropy of the plasma by measuring the rapidity de- 
pendence of photon production at various energies. 



PHOTON RATE - HARD CONTRIBUTION 

There are two hard contributions coming from the an- 
nihilation and Compton scattering diagrams. The rate of 
photon production from the medium quark annihilation 
diagrams can be expressed as 
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where we have assumed massless quarks, the Mandel- 
stam variables are defined hy t = {Ki — Q)^ and u = 
{K2 — Q)^, — 5/9 comes from SU(3) Casimirs, and 
= J (Pk/{2TT)^. To make the distinction between 
four-vectors and three-vectors clear we use upper-case 
letters for four- vectors. Lower-cased letters will repre- 
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sent the length of a three- vector while a bold-faced lower- 
cased letter stands for a three-vector. Here fq^g are the 
medium quark and gluon distribution functions and we 
have assumed that the distribution functions for quarks 
and anti-quarks are the same, fg = fq. 

In this work we will also assume that the anisotropic 
distribution functions are obtained by squeezing or 
stretching an arbitrary isotropic distribution along one 
direction in momentum space, i.e., /i(k, ^,phard) = 

/i.iso ("v/k^ +C(k- n)^,Phard) , where i = {q,g}, Phard is 

a hard momentum scale which appears in the distribu- 
tion functions, n is the direction of the anisotropy, and 
^ > — 1 is a parameter reflecting the strength and type of 
the anisotropy [4]. In the following we will suppress the 
explicit dependence on ^ and Phard- 

Anticipating the infrared divergence associated with 
this graph we begin by first changing variables in the 
first integration to P = Ki — Q and introduce an in- 
frared cutoff p* on the integration over the exchanged 
three- momentum p [24] . We then choose spherical coor- 
dinates with the anisotropy vector fi defining the z-axis. 
Exploiting the azimuthal symmetry about the 2;-axis we 
also choose q to lie in the x — z plane. Using the delta 
function to perform four out of the nine integrations and 
expanding out the phase-space integrals explicitly gives 
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with t = u)'^ — p'^ , u ^ {k - q)^ — {k — q)^, and w = 
|p + q| — q. The azimuthal angles (pi are defined through 
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2pk sin 9p sin 6k 
and X > is given by 

X = 4p^A;^sin^^fesin^^p- 
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The rate of photon production from the Compton scat- 
tering diagrams can be expressed as 
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where we have again assumed massless quarks, fq = fq, 
and the Mandelstam variables are defined by s = {Ki + 
and t = {Ki—Q)^. Performing a change of variables 
to P = Ki — Q and proceeding as with the annihilation 
contribution gives 
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where here t = ui^ - p^ , s = {oj + q + k)"^ - {p + q + k)^, 
and recalling lu = \p + q\ — q. 

The total photon production from processes which 
have a hard momentum exchange is given by the sum 
of Eqs. (2) and (6). 
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To evaluate the five-dimensional integrals in Eqs. (2) and 
(6) we use Monte-Carlo integration. The total hard con- 
tribution (7) has a logarithmic infrared divergence as 
p* — !■ 0. This logarithmic infrared divergence will be 
cancelled by a corresponding ultraviolet divergence in the 
soft contribution. 



PHOTON RATE - SOFT CONTRIBUTION 

We now turn to the calculation of the previously ex- 
cluded contribution involving soft momentum exchange. 
We adopt the Keldysh formulation of quantum field the- 
ory, which is appropriate for systems away from equilib- 
rium [25, 26]. In this formalism both propagators and 
self energies have a 2 x 2 matrix structure. The compo- 
nents (12) and (21) of the self-energy matrices are related 
to the emission and absorption probability of the particle 
species under consideration [27-29]. Due to the almost 
complete lack of photon absorbing back reactions, the 
rate of photon emission can be expressed as [30] 
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from the trace of the (12)-element II12 of the photon- 
polarization tensor. 

We evaluate Tli2'^ using the hard-loop (HL) resummed 

fermion propagator derived in Ref. [31]. In that Refer- 
ence we showed that the needed off-diagonal components 
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of the fermion self-energy can be expressed in terms of 
the retarded self-energy 
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where 

/(k)^2(/,(k)+/,-(k))+4/,(k). (10) 
Taking the HL limit where appropriate we find that 

»ni2{:(g) = -e^e^iV.^AM f g,A'^(p) , (11) 
where 
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The tensor A is defined through 
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where the star indicates complex conjugation. To com- 
pute the soft photon rate (8) wc must specify the 
anisotropic distribution functions we use and then nu- 
merically evaluate Eqs. (9) and (11) with an ultraviolet 
cutoff p* placed on the length of three-momentum. 



RESULTS 

The total photon rate is given by adding the hard 

(7) and soft (8) rates. Here wc specialize to the case 
^ > which is the relevant one for times r > 0.2 fm/c. 
For the arbitrary isotropic distribution functions, /i^iso, 
which are contracted along the z-axis, wc choose Bose- 
Einstein and Fermi-Dirac distributions in the case of glu- 
ons and quark/anti-quarks, respectively. Using these dis- 
tribution functions wc reproduce the equilibrium results 
which have been obtained previously [18, 19] by taking 
the limit ^ ^ 0. For ^ ^ we find that the soft /hard 
logarithms of p* cancel for all photon angles/energies. 

In Fig. 1 we plot the dependence of the total 
anisotropic photon production rate as a function of pho- 
ton energy, E, for three different photon propagation an- 
gles assuming ^ = 10 and as = 0.32. This Figure shows 
that there is a clear dependence of the spectrum on pho- 
ton angle with the difference increasing as the energy of 
the photon increases. The shaded bands indicate our es- 
timated theoretical uncertainty which is determined by 
varying the hard-soft separation scale p* by a factor of 
two around its central value which is parametrically p* ~ 
\/9Pha.Td- We have scaled everywhere by the arbitrary 
hard-momentum scale phard which appears in the quark 
and gluon distribution functions. This scale will be time 




FIG. 1: Photon rate for ^ = 10 and = 0.32 as a function of 
energy for three different photon angles, 9q — {0, 7r/4, 7r/2}, 
corresponding to rapidities y = {oo, 0.88, 0}, respectively. 
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FIG. 2: Photon rate for ^ = {0, 1, 10} as a function of the 
photon angle, 6q, for S/phard = 5 and as = 0.32. 



dependent with its value set to the nuclear saturation 
scale, Qs, at early times and to the plasma temperature 
at late times. For RHIC ^ 1.4—2 GeV and expected 
initial plasma temperatures are To ^ 300 — 400 MeV. 

To compare the angular dependence at different ^ in 
Fig. 2 we show the dependence of the photon produc- 
tion rate on angle at a fixed photon energy i^/phard = 5 
and as = 0.32. As can be seen from Fig. 2 the dif- 
ference between the forward and transverse production 
rates increases as ^ increases. To summarize the ef- 
fect we define the photon anisotropy parameter, TZ^ = 
(dR/d^qle^^o) / (dR/d^q\g^^T^/2) ■ This ratio is one at all 
energies if the plasma is isotropic. For anisotropic plas- 
mas it increases as the anisotropy of the system and/or 
the energy of the photon increases. Due to the limited 
experimental rapidity acceptance one could define this 
ratio at a lower angle, e.g. 7r/4. However, the ability to 
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TABLE I: Photon anisotropy parameter TZ-^ for different pho- 
ton energies and ^, assuming as = 0.32. 



resolve anisotropies increases as the sensitivity to forward 
angles increases so it is best to compare the most forward 
photons possible with transverse photon emission. 



CONCLUSIONS AND DISCUSSION 

The results presented in the previous section demon- 
strate that the differential high-energy medium pho- 
ton production rate is sensitive to quark-gluon plasma 
momentum-space anisotropies with the sensitivity in- 
creasing with photon energy. Turning this into the final 
expected experimental yields as a function of photon en- 
ergy and angle will require folding this rate together with 
a phcnomenological model of the time dependence of ^ 
and Phard- Current models of plasma evolution do not 
take into account the evolution of the momentum-space 
anisotropy of the parton distribution fimctions. They 
therefore implicitly assume an isotropic thermal plasma 
at all times so it will be necessary to extend these models 
to include the time dependence of ,j and Phard- However, 
since high-energy medium photon production is domi- 
nated by early times these photons will be mostly sensi- 
tive to the first 1-2 fm/c of plasma evolution where simple 
models will hopefully be efficacious. 

One theoretical shortcoming of our calculation is that 
it does not include the bremsstrahlung contribution 
which in an equilibrium plasma also contributes at lead- 
ing order in the coupling constant due to enhancements 
from coUinear photon radiation [20 23]. This contribu- 
tion has been omitted because it is currently not known 
how to resolve the problem that the presence of unstable 
modes causes unregulated singularities to appear in ma- 
trix elements which involve soft gluon exchange. Recent 
work [32] has suggested that these singularities could be 
shielded by next-to-leading order corrections to the glu- 
onic polarization tensor, however, the detailed evalua- 
tion of these NLO corrections has not been performed 
to date. Absent an explicit calculation our naive expec- 
tation would be that these NLO corrections would yield 
similar angular dependence since they are also peaked at 
nearly coUinear angles. 
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